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The anharmonic electron-phonon problem is solved in the 
infinite-dimensional limit using quantum Monte Carlo simula- 
tion. Charge-density-wave order is seen to remain at half fill- 
ing even though the anharmonicity removes the particle-hole 
symmetry (and hence the nesting instability) of the model. 
Superconductivity is strongly favored away from half filling 
(relative to the charge-density-wave order) but the anhar- 
monicity does not enhance transition temperatures over the 
maximal values found in the harmonic limit. 



Introduction Lattice anharmonicity is responsible for 
many different properties of solids. For example, the ex- 
pansion of a crystal upon heating is due entirely to an- 
harmonic effects — harmonic crystals do not change their 
volume upon heating. In addition, the phonon contri- 
bution to thermal resistivity is due, in some measure, 
to anharmonic effects. Furthermore, lattice modes in a 
real crystal are never completely harmonic. Anharmonic 
terms always enter to higher-order when expanding in 
the amplitude of the lattice vibrations. 

But the effects of anharmonicity on the electron- 
phonon problem, and more specifically on superconduc- 
tivity, are relatively unknown. This is surprising, because 
the electron-phonon problem is well studied, especially in 
the weak coupling limit. Here the theory of electrons in- 
teracting with harmonic phonons, initially pioneered by 
Migdal J!]] and Eliashberg has led to an accurate de- 
scription of superconductivity in most low-temperature 
phonon-mediated superconductors. It has been found, 
that in the regime where the electron-phonon coupling is 
weak so that the Migdal-Eliashberg theory still applies 
to the anharmonic case, the superconducting transition 
temperature is not enhanced by the presence of lattice 
anharmonicity |3|. The situation as the coupling is in- 
creased is less clear. Recently, a new mechanism for 
superconductivity, driven by kinetic-energy effects, has 
been hypothesized in anharmonic crystals that display 
stronger electron-phonon coupling ^] . In this new mech- 
anism, paired holes, in a nearly filled band, are more 



mobile than single holes, driving the superconductivity. 

An exact treatment of lattice anharmonicity is, how- 
ever, difficult from a theoretical point of view because an 
anharmonic "perturbation" is never small — the phonon 
wavefunctions are always dominated by the anharmonic 
terms in the potential as the phonon coordinate becomes 
large. All that is known rigorously about the anharmonic 
electron-phonon problem is that the ground state must 
contain a spin singlet Q for even numbers of electrons 
on a finite lattice. 

Much progress can be made in the limit of large spa- 
tial dimension |(| where the lattice many-body problem 
can be mapped onto a self-consistently embedded impu- 
rity problem that is solved via a quantum Monte Carlo 
(QMC) simulation. 

Model The simplest electron-phonon model that in- 
cludes anharmonic effects is the anharmonic Holstein 
model Q in which the conduction electrons interact with 
local phonon modes: 
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where cj CT (cj CT ) creates (destroys) an electron at site i 
with spin a, n- la = cl a Ci a is the electron number opera- 
tor, and Xi (pi) is the phonon coordinate (momentum) at 
site i. The hopping of electrons is restricted to nearest- 
neighbor lattice sites i and j and is governed by the hop- 
ping matrix element Uj (tij is a Hermitian matrix). 

The local phonon has a mass M and a frequency £1 
associated with it. The anharmonic contribution to the 
phonon potential energy is chosen to be a quartic in the 
phonon coordinate with a strength a an h- The deforma- 
tion potential (electron-phonon interaction strength) is 
parameterized by an energy per unit length and is de- 
noted g. A useful combination of fundamental param- 
eters is the bipolaron binding energy (in the harmonic 
limit) U := — g 2 /Alii 2 which determines the energy scale 
for the effective electron-electron interaction mediated by 
the phonons. The chemical potential is \x. 

The hopping ty is isotropic with strength t =: t* /2y/d 
where t* = 1 to define the energy scale and d is the 
dimensionality. The mass is then set equal to one (M = 
1) leaving U, a an h, £1, and the electron filling as free 
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parameters. We set ft = 0.5t* which is approximately an 
order of magnitude smaller than the effective bandwidth 
(W w it*), and we set U = -it* (g = 1) which is well 
into the strongly coupled regime in the harmonic limit. 

Results The original Holstein Hamiltonian corre- 
sponds to the case a an h = 0. Both the harmonic Holstein 
model and the harmonic Holstcin-Hubbard model have 
been solved exactly in the limit of infinite dimensions via 
quantum Monte Carlo simulation dH|. These models 
display charge-density-wave (CDW) order near half fill- 
ing and superconductivity (SC) away from half-filling. As 
the phonon frequency is increased, the SC is favored rela- 
tive to the CDW order. However, in the strong-coupling 
limit, CDW order is favored over SC because of the band- 
narrowing effect of the bipolaron. The quantum Monte 
Carlo simulations also found that the effective phonon 
potential (determined after integrating out the effects of 
the electrons) generically acquires a double- well structure 
signifying the high-temperature formation of a bipolaron. 

An initial analysis of the anharmonic model can be 
made by employing the Born-Oppenheimer approxima- 
tion in the strong-coupling limit (Uj = 0) B : the phonon 
frequency is assumed to be smaller than any of the other 
energy scales so the phonons can be approximated by 
static lattice distortions corresponding to the minimum 
of the local phonon potential energy. Since the phonons 
couple linearly to the electronic charge, the equilibrium 
phonon coordinate varies when there are zero, one, or 
two electrons on a site. The origin xq = is chosen to 
correspond to the case with no electrons on a site. Then 
x\ and X2 denote the equilibrium coordinates with one or 
two electrons on a lattice site. In the harmonic case the 
relative distances Xi — x$ and X2 — x\ are identical, which 
is a requirement for particle-hole symmetry. When a lat- 
tice anharmonicity is turned on, the equilibrium phonon 
coordinates with one and two electrons on a lattice site 
all move toward the origin, but the relative distances are 
no longer symmetric, rather the distance xi — x\ becomes 
significantly smaller than x± — xq, as can be seen in Fig- 
ure 1(a). 

The lattice anharmonicity and subsequent asymmetric 
shift of the equilibrium phonon coordinate have two main 
effects: (I) the model loses particle-hole symmetry which 
may allow a new type of kinetic-energy driven super- 
conductivity to emerge j| and (2) the effective electron- 
electron attraction is dramatically reduced as can be seen 
by a plot of the bipolaron binding energy in Figure 1(b). 
Thus the lattice anharmonicity generates an effective re- 
tarded repulsive interaction between the electrons. It 
is possible that this effective repulsive interaction be- 
tween the electrons is sufficiently strong to destroy the 
kinetic-energy driven superconductivity. The lattice an- 
harmonicity also breaks particle-hole symmetry, remov- 
ing the nesting instability of the CDW at half filling and 
weak coupling. Hence, in the weak-coupling limit, one 
expects the lattice anharmonicity to favor SC relative to 



CDW order, possibly completely suppressing CDW or- 
der, although it is also likely that anharmonic effects will 
reduce all transition temperatures. What is surprising 
is that a lattice anharmonicity as small as a an h — 0.01 
(cf. Fig. 1(b)) can have such a large effect on the electron- 
phonon problem in the strongly coupled regime studied 
here. 

The anharmonic Holstein model is solved exactly in 
the infinite-dimensional limit via QMC simulation. The 
algorithm is straightforward, and has been described nu- 
merous times before |n]||[)| , so it will not be repeated 
here. The QMC algorithm determines the single and 
two-particle Green functions at each of the Matsubara 
frequencies, which allows one to construct directly the 
momentum-dependent susceptibility to either CDW or 
SC order of the infinite-dimensional lattice. Divergence 
of either of these susceptibilities signals a transition to 
an ordered state. 

We find that at half filling the system orders in a com- 
mensurate CDW phase. The results of the QMC simu- 
lation are presented in Figure 2. The solid dots repre- 
sent the QMC transition temperatures, and the dashed 
line is a guide to the eye. The solid line is the result 
of a strong-coupling perturbation theory expansion |ll| ] 
(fourth-order perturbation theory in t*) that predicts 
that the system evolves from CDW order (solid line) to 
SC order (dotted line). We find no evidence for such 
an evolution in the QMC solution. Rather, the system 
remains in the CDW-ordered state even though there is 
no particle-hole symmetry remaining in the model. It is 
possible that as the strength of the anharmonicity is in- 
creased further that the SC order may overtake the CDW 
order, but it would be difficult to see that behavior with 
the QMC simulations because the transition tempera- 
tures become too small to be reliably determined. Note 
that the strong-coupling perturbation theory is accurate 
in predicting the peak position and peak height for the 
maximal CDW transition temperature. This is similar 
to what happened in the harmonic case as a function of 
the coupling strength. 

The initial rise in T c occurs because of the reduction in 
the bipolaron binding energy E^p as the strength of the 
anharmonicity increases. In the strong-coupling regime, 
the transition temperature is inversely proportional to 
Ebip (T c (x t* 2 /Ei,ip). As the bipolaron unbinds, the 
system evolves from preformed-pair physics to BCS-type 
physics where pair formation and condensation both oc- 
cur at T c . This then causes the transition temperature 
to drop as a an h is increased further. 

This physics can be dramatically illustrated by ex- 
tracting the effective phonon potential from the QMC 
simulation M . The probability distribution for the local 
phonon coordinate x is measured for each of the imag- 
inary time slices. This probability distribution is then 
equated to an activated form for an effective potential 
P{x) = exp[— V e f f [x)/T] . The effective phonon potential 
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includes the effects of the interaction of the electrons with 
the phonons. In the harmonic limit, this effective poten- 
tial was seen to evolve from a harmonic form [x — x\) 2 
to a symmetric double- well form as the coupling strength 
increased 11 . The double well signified bipolaron forma- 
tion, since the number of electrons sitting at the lattice 
site was either zero or two for each of the minima (cen- 
tered around xg and X2 )• The results for the effective 
phonon potential in the anharmonic model at half fill- 
ing and T = 0.125i* are presented in Figure 3. The 
effective potential is plotted versus a rescaled phonon co- 
ordinate for different values of a a nh- The zeroes of the 
potentials are shifted vertically to separate the different 
plots. Note that the addition of anharmonicity squeezes 
the phonon coordinate toward the origin, and makes the 
effective phonon potential asymmetric (corresponding to 
the loss of particle-hole symmetry). The barrier height 
decreases as a an h increases, as required by the reduction 
of the bipolaron binding energy. Note that the maximum 
CDW T c occurs when the barrier height is on the order 
of T c (a anh « 0.0075). 

In order to find SC solutions, we must dope the system 
off half filling. The phase diagram for three different val- 
ues of a an h is shown in Figure 4. In Figure 4(a), a small 
value of a an h = 0.005 is chosen. This result is close to the 
harmonic result. Nevertheless, even at this small value of 
oianhi the phase diagram is not symmetric with respect 
to reflection about half filling, signifying the particle-hole 
asymmetry. We compare to the strong-coupling pertur- 
bation theory, which predicts CDW order to survive until 
doped well off of half filling. This approximation is, by 
construction, particle-hole symmetric. Note that the case 
of less than half filling is approximated much better by 
the perturbative calculation, than the greater than half- 
filling case. 

Figure 4(b), has a an h = 0.02 which lies at a point be- 
yond where the CDW T c is maximal. The perturbative 
approximation is less accurate here, but displays the cor- 
rect order-of-magnitude estimate for the phase diagram. 
We see how the anharmonicity favors the superconduct- 
ing solutions, pushing the CDW-SC phase boundary to- 
ward half filling. SC is favored because the squeezing of 
the phonon coordinate toward the origin enhances the 
Frank-Condon overlap factors associated with electronic 
hopping and increases the electronic bandwidth. 

Figure 4(c), takes the large value a an h = 0.1. The 
strong-coupling theory fails in this regime, and so it has 
not been included in the plot. We see once again how the 
increase in a an h now suppresses the transition tempera- 
tures, and how it moves the CDW-SC boundary closer to 
half filling. There is no evidence for the new SC mecha- 
nism. 

We found no evidence for incommensurate CDW order, 
either at half filling or away from half filling. The CDW 
susceptibility still had it's maximum at the (tt, tt,tt, . . .) 
point when the SC susceptibility diverged at a higher T c . 



Anharmonicity does not favor incommensurate order in 
this regime of parameter space. 

The effect of anharmonicity on the electron-phonon 
problem has some aspects that are similar to, and some 
aspects that are different from the effect of a direct 
Coulomb repulsion |9). In both cases the superconduct- 
ing order is favored by the additional interaction rela- 
tive to the charge-density-wave order, in the sense that 
the SC-CDW phase boundary moves toward half fill- 
ing. However, the anharmonicity breaks the particle- 
hole symmetry, while the Coulomb repulsion does not. 
Furthermore, the CDW transition temperature is moder- 
ately suppressed even for large values of the anharmonic- 
ity. Whereas, if the Coulomb repulsion is increased to 
large enough values, the CDW T c will plummet to zero. 

Conclusion We have analyzed the effect of anhar- 
monicity on the electron-phonon problem in the regime 
where the harmonic model has formed bipolarons. We 
find that the anharmonicity breaks particle-hole symme- 
try and produces phase diagrams that are not symmetric 
about half filling. The anharmonicity acts in some re- 
spects as a retarded "Coulomb repulsion" , since it greatly 
favors SC solutions relative to CDW ones. However, in 
the infinite-dimensional limit, we did not find any en- 
hancement to T c 's for the anharmonic model relative to 
the maximal values found in the harmonic case. Fur- 
thermore, we found no evidence for any kind of novel 
superconducting mechanism driven by kinetic-energy ef- 
fects. We also studied the effective phonon potential 
which includes the effects of the interaction of the phonon 
with the electrons. We discovered that the anharmonicity 
tends to squeeze the phonon coordinates closer to the ori- 
gin, and rapidly suppresses the bipolaron formation. We 
also verified that the strong-coupling perturbation the- 
ory accurately describes the anharmonic Holstein model 
in the regime where the system still has preformed pairs. 
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FIG. 1. (a) Phonon coordinates at the minima of the an- 
harmonic phonon potential with zero (xo), one (xi), and two 
(22) electrons on a site as a function of a an h in the atomic 
limit (tij = 0) . (b) Reduction of the bipolaron binding energy 
as a function of a an h in the atomic limit (tij = 0) . 
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FIG. 2. Charge-density-wave transition temperature at 
half filling as a function of ct anh . The QMC results (dots-no 
SC transitions were found at half filling with the QMC) 
are compared to the strong-coupling perturbation theory for 
CDW order (solid line) and SC order (dotted line). The 
dashed line through the QMC points is a guide to the eye. 



FIG. 3. Effective phonon potential for the anharmonic Hol- 
stein model at half filling and T = 0.125i*. The labels on the 
curves are the values of a anh . The potentials are shifted by 
adding a constant to them to separate the curves. 



FIG. 4. Phase diagrams for the anharmonic Holstein model \ 

at three different values of ctanh'- (a) 0.005; (b) 0.02; and (c) * x 

0.1. The solid dots are CDW phases, the open triangles are SC it 

phases. The solid lines are a strong-coupling approximation > 
to the CDW T c , and the dotted lines are a strong-coupling 
approximation to the SC T c . 




0.05 



CDW (approx) — 
SC (approx) - ■ 
CDW (qmc) . 
SC (qmc) A 



_l_ 



_l_ 



_l_ 



0.05 0.1 0.15 0.2 

«anh 

Freericks, Jarrell, and Mahan, Phys. Rev. Lett., Figure 2 




x*=xM[r/2g 

Freericks, Jarrell, and Mahan, Phys. Rev. Lett., Figure 3 



4 




CDW (approx) - 
CDW (qmc) < 



(o) 



w 02 

o 

° 0.15 



V 0.05 



^ 0.2 
o 

II °- 15 
I 0.1 

V 0.05 



Electron concentration p 



I 

SC (approx) - - - 


CDW (approx) (b) 


- SC (qmc)A 


CDW (qmc) • 










l aa a * 


• 

A A , 


1 


, A A 



0.5 1 1.5 

Electron concentration p e 



1 

SC (approx) - - - 


CDW (approx) (c) 


- SC (qmc)A 


CDW (qmc) • 








AA "V A i A A 


) 0.5 


1 1.5 ; 



Electron concentration p e 
Freericks, Jarrell, and Mahan, Phys. Rev. Lett., Figure 4 



